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Abstract 

Electron-deuteron elastic scattering data ( A(Q 2 ) and B[Q 2 ) struc¬ 
ture functions and polarization observables t 2 o, t 2 ± and £ 22 ) are fit 
with a model that respects asymptotic properties of pQCD at high 
momentum transfer. The data analysis shows that pQCD starts from 
Qqcd = 3-5 (GeV/c) 2 . Predictions for the magnetic structure func¬ 
tion B(Q 2 ) and the polarization observables at high momentum trans¬ 
fer are given. 

PACS number(s): 13.40.Fn, 12.38.Bx, 25.30.Bf 


1 Introduction 

In recent few years new data from TJINAF on the ed elastic scattering were 
reported. They include the electric structure function, A(Q 2 ), measured with 
high precision up to Q 2 = 6 (GeV/c) 2 [I]. |2J and measurements of tensor 
polarization observables, t 2 0, £21 and t 2 2, up to Q 2 = 1.7 (GeV/c) 2 ||. 

This data, together with data on the magnetic structure function, B(Q 2 ) 
|| , restrict the deuteron structure at scales where quark-gluon degrees of 
freedom are expected to become defrozen. For example, according to opti¬ 
mistic estimations pQCD should start from Q 2 of order of few (GeV/c) 2 J4|. 
It is nice that this prediction was confirmed by analysis of TJINAF data on 
A{Q 2 ) at Q 2 > 2 (GeV/c) 2 [|. 

For further conclusions one also should consider the spin structure of the 
deuteron from pQCD point of view. However data on polarization observ¬ 
ables, as well as on B(Q 2 ), correspond to Q 2 < 2 (GeV/c) 2 , which is not 
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enough for pQCD. This is a typical intermediate region between nucleon- 
meson and quark-gluon pictures, where isobar configurations, meson ex¬ 
change currents and constituent quark degrees of freedom are all important 

i- 

The purpose of this work is to investigate phenomenologically a smooth 
connection between nucleon-meson and pQCD regions and make predictions 
for B(Q 2 ) and the polarization observables at higher Q 2 , where pQCD should 
work. A parameterization which connects these two regions was proposed 
earlier by one of the authors (A.P. K.) and A.I. Syamtomov |7j. It assumes 
power fall off of hclicity spin amplitudes at asymptotically high Q 2 coming 
from quark counting rules. A new analysis of the parameterization 0 which 
includes the recent TJINAF data was provided in Now we study loga¬ 
rithmic corrections to the power behavior. Such corrections are shown to be 
important for the structure function A(Q 2 ) at the highest region of TJINAF 
energy @. 

The paper is organized as follows. In sect. 2 we discuss the general struc¬ 
ture of the helicity amplitudes for the elastic ed scattering in the light cone 
frame (LCF) and pQCD predictions for the hclicity amplitudes at high Q 2 . 
Parameterization of the hclicity amplitudes which smoothly connects regions 
of low and high Q 2 is given in sect. 3. Then, in sect. 4, the data base and 
fitting procedure are summarized. Discussions and summary are given in 
sect. 5. 

2 Helicity amplitudes and the deuteron form 
factors 

2.1 Helicity amplitudes in LCF 

The main object of our analysis is the hclicity amplitudes of the 7 * + d —> d 
transition 

= (P, A>, (1) 

where p and p' are momenta and A and X 1 are hclicities of the deuteron in 
the initial and final states, respectively. 

Due to gauge invariance, covariance and discrete symmetries only three 
of the 36 hclicity amplitudes (|l|) are independent and one can choose dif¬ 
ferent sets of independent helicity amplitudes. Direct calculations, however, 
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demonstrate that it is not so in dynamics at LCF This phenomena was 
shown to come from the incompatibility of transformation properties of ap¬ 
proximate current and a deuteron wave function used in practical calculations 


10 , 11]. As a result a non-physical dependence on orientation of light-front 


plane appears. Thus the choice of the independent amplitudes becomes of 
great importance in pQCD calculations where LCF is often used. 

First let us define LCF as a limiting reference system where the ^-projection 
of the incoming and outgoing deuteron is close to infinity. In LCF the mo¬ 
menta of the incoming and outgoing deuteron are given as follows 


-- - 


( 2 ) 


with 


p + = p° + p 3 = 2V, p' + = p° + p' 3 = 2V, VT>M 2 +p 2 ± (3) 

(M is the deuteron mass). The momentum of the virtual photon is given by 

( 4 ) 


= (0,—2ph,0), £l = -( -Q,0 

and the polarization vectors for the deuteron in the initial and final states, 
respectively, read 
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Here we put p± = Using the standard expression for the e.m. current 

matrix element 


Jyx= ~ {Gi W 2 )£’(A\ p')e(\, p)(p + p')"+ 

+ G 2 p) (e"(\',p')q) (e(A,p)j)]- 

_ G , w . )(p+j0 


( 6 ) 
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one gets the following expressions for the current plus-component 

J+ = p + {2(l-2/ / )G 1 + 4r / G 2 -4p 2 G 3 }, (7) 

J [o = P + 12 \J~2rjG i — x /2j]G 2j + 2a^/2^7^7^3 j - , (8) 

Jt-i = ~P + {2 V G 3 }, (9) 

J £L = P + {2G! + 277G3} (10) 

where J A J“ A , = J AA , + J AA ,. It is easy to show that they satisfy the so-called 
angular condition 

(1 + 2 7j)J+ + J+ ! - 2^7+ - J+ = 0 (11) 


and thus there are only three independent helicity amplitudes between the 

J\ii J\-\ 1 -7io an( l 4o B P~^1 • 

Alternatively the angular condition (|TT|) teaches us that even at pQCD 
extreme there appears (through dimensionless ratio p = ^ 2 ) an additional 
scale parameter 4M 2 , apart from the pQCD parameter Aq CD . 

2.2 The deuteron form factors 

The charge, Gc(Q 2 ), magnetic, Gm(Q 2 ), and quadruple, Gq(Q 2 ), form fac¬ 
tors are connected with the form factors G\(Q 2 ), G 2 (Q 2 ) and G 3 (Q 2 ) as 
follows 


Gq — G\ — G 2 + (1 + p)G 3 , 

Gc = G\ + -t/Gq, 

Gm = G 2 . ( 12 ) 

Using (0)-(|lO;) one expresses Gc{Q 2 ), Gm{Q 2 ) and Gq(Q 2 ) in terms of any 
three helicity amplitudes J xx ,, for example 

G c = 

Gm = 

Gq = 


2V{2p + 1) 

1 

2V(2rj + 1) _ 
1 

2V(2 V + 1) 
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6 ' yoo + 3 V 2 Jlo+ 3 '- 1 


( 13 ) 
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In turn, the A(Q 2 ) and B(Q 2 ) structure functions and the t 2 i( 6 ,Q 2 ) polar¬ 
izations read 
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77 (77 + 77 sin - G m Gq , 
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where S = d 2 + S 2 tg 2 §. 


(14) 

(15) 

(16) 

(17) 

(18) 


2.3 pQCD predictions for the helicity amplitudes 


From pQCD arguments one gets very simple rules to determine the power 
behavior of the helicity amplitudes J^, J ± 0 and J+_ x [|13|, [L4j]. For example, 
it follows that the amplitude is a leading amplitude with an asymptotic 
fall off 

' A f -) 10 w 

up to logarithmic corrections. It was also argued that in LCF the helicity 
flip amplitudes Jf Q and J [Li are suppressed as [[15] 
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( 20 ) 


Similar considerations give that 
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which agrees with the angular condition (|TTD , but only at extremely high Q 2 , 
Q 2 > AM 2 . 
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In our analysis, following arguments of Ref. [|T3|, |T(|, we consider the 
set J^, Jf 0 and J+_ ± as main amplitudes, where the behavior (|T9|) , (| 20 |) is 
regulated in the intermediate region only by the Aqcd- In turn, the amplitude 
Jf 1 ) must be determined from the angular condition (|IT|) and thus it depends 
on the two scale parameters, Aqcd and AM 2 . 


3 Parameterization of helicity transition am¬ 
plitudes 


Following the idea of reduced nuclear amplitudes in QCD [|17]], we define the 
reduced helicity transition amplitudes goo, go+ and g + - as follows: 


A p JtAQ 2 ) = a 2 (^Q 2 '^9,AQ 2 ), 


( 22 ) 


where G(Q 2 ) is a three-quark-cluster (nucleon) form factor. For the G(Q 2 ) 

I -2 


we assume dipole behavior G(Q 2 ) — 1 + yr , but with the parameter /jl 2 

is different from that for a free nucleon 0.71 (GeV/c) 2 . 

We consider separately two kinematical regions, large Q 2 region, Q 2 > 
Qqcdj an d low Q 2 region, Q 2 < Q q CD . The parameter Qq CD is expected to 
be of order of few (GeV/c) 2 . Its exact value will be determined from a fit to 
experimental data. 

pQCD predicts that at asymptoticaly high Q 2 the reduced transition 
amplitudes behave as follows 

n s£ y) = 


(asy) 

%o = 

Ah 

Q 2 

(asy) 

n 3 

9+- = 

Q A 


In (|23|) the factor 4>{Q 2 ) takes into account logarithmic corrections 

[a,(Q 2 )] 5 (lgQ 2 /A|oDF 


<P(Q ) = 


K(Q 2 /4)] 1 [lg« 2 /(4A|cD)K 


(23) 


(24) 


and and 7 ^ are leadings anomalous dimensions for the deuteron and the 
nucleon, respectively, 

^ 6 Cf n Cp 


7 = 


5 P' 1 2/3 ’ 


(25) 
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where Cp = [n 2 — l)/( 2 n c ), /5 = 11 — |n/, n c = 3 is number of quark colors 
and nf — 2 is the number of flavors, a s (Q 2 ) = j3 ln g 2 / A a is the running 

quark-gluon coupling Q. 

At Q 2 < Qqcd the following parameterization for the reduced amplitudes 
is assumed 


9 oo 
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where a 2 


[ck 0 + (n — l)m 0 ] 2 . One imposes the constrains (see, e.q., [0): 
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on the coefficients a n , b n and c n to demand form factor normalization at 
Q 2 = 0. In (g7|) /j,d = 0.857406 M/m p is the deuteron magnetic moment in 
“deuteron magnetons” and Qd = 25.84 is the deuteron quadrupolc momen¬ 
tum in M 2 je. 

The coefficients Ni, N 2 and 7V 3 appearing in (|23|) cannot be calculated 
from pQCD and we determine them from the smooth connection of the two 
parametrizations at the point Q 2 = Qqqd 
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4 Data base and fitting parameters 


In out fit we use the following data: for A(Q 2 ) from |], |2], [18], [A]. |2T], [22], ^3|, [24|. 
. for B(Q 2 ) from [g [A], P, |6|, |7| , and for t 20 from @, @ p. |34 


and for t 2 2 , t 2 i from ||. [Il|]. But it must be noted that due to large errors 
t 22 , t 2 1 data are practically not informative for the fit. Data for A from 
27] are in sharp contradiction with all the wold data set and therefore were 
omitted from the fit. We also do not include in the data base results from |35 


because they have large errors and practically do not change the obtained 
results. Data |23| are somewhat lower that the world data, but due to large 


their uncertainties the fit is insensitive to this data. Data for the B structure 
function were normalized, if necessary, to our convention of the magnetic 
form factor. 

Qqcd was considered as a parameter of the model, but the QCD cutoff 
parameter Aq CD was fixed at 200 MeV. In (|27l) we chose N — 4, so that the 
the model has 10 independent fit parameters. The fit parameters summarized 
in Table 1 were obtained with x 2 — 399 for 200 data points. 


5 Discussion and summary 

Figs. 1 and 2 display comparison of the model with A(Q 2 ) data. To show 
how far experimental data are from pQCD results in preasymptotic region 
we continue asymptotic behavior given by (|23|) to region lower than Q q CD 
(dotted line on Fig. 2 and Figs. 3-7). One concludes that for A(Q 2 ) pQCD 
works from Q 2 -region between 1 and 2 (GeV/c) 2 . Comparison with data for 
B(Q 2 ) and the polarization observables are given on Figs. 3 — 5. One sees 
that for the magnetic form factor pQCD should start from Q 2 between 2 
and 3 GeV/c, but for the polarization observable it should start somewhat 
earlier, near 2 (GeV/c) 2 . 

Fig. 7 shows results for for the charge and quadrupole form factors, 
G C (Q 2 ) and G q (Q 2 ). 

In summary, we give a parameterization of the deuteron form factors up 
to Q 2 = 6 (GeV/c) 2 . Asymptotic behavior of the form factors is dictated 
by quark counting rules and pQCD helicity rules and therefore one may 
hope that this parameterization can be extrapolated for higher transferred 
momentum. For example, the model predicts behavior of the magnetic struc¬ 
ture function, B(Q 2 ), at Q 2 > 2.5 (GeV/c) 2 which can be studied in future 





























































experiments. 
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Table 1: Parameters of the model. 


OL 0 

0.2635 GeV 

m 0 

0.36864 GeV 

ai 

0.5373059 • 10" 1 GeV~ 2 

0,2 

0.2836982 GeV“ 2 

03 

from eqs. (28) 

04 

from eqs. (27) 

h 

-9.898079- 10” 1 

h 

1.6950028 

h 

from eqs. (28) 

h 

from eqs. (27) 

Cl 

-0.1074219 

c 2 

-0.7075449 • 10" 1 

c 3 

from eqs. (28) 

c 4 

from eqs. (27) 

lllFiBfi 

3.50120 GeV 2 

n' 2 

0.50959 GeV 2 
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Figure captions 


Figure 1. Comparison of the model fit with data for A(Q 2 ) at Q 2 < 
1 (GeV/c) 2 . 

Figure 2. Comparison of the model fit (solid line) with data for A(Q 2 ) at 
Q 2 > 1 (GeV/c) 2 . The dotted line is asymptotic behavior given by d23|) 
extrapolated to lower transfer momentum. 

Figure 3. Comparison of the model fit with data for B(Q 2 ). The lines are 
the same as at Fig. 2. 

Figure 4. Comparison of the model fit with data for f 2 o- The lines are the 
same as at Fig. 2. 

Figure 5. Comparison of the model fit with data for f 2 i • The lines are the 
same as at Fig. 2. 

Figure 6. Comparison of the model fit with data for t 2 2 - The lines are the 
same as at Fig. 2. 

Figure 7. Comparison of the model fit with data for |G C (Q 2 )| (upper figure) 
and G q (Q 2 ) (lower figure). Data are from ||. For the last two points the 
two solutions (filled and open points) are shown, see [§. The lines are the 
same as at Fig. 2. 
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